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a generalization of the Chirikov- Taylor model is introduced to study the dynamics 



< 



of a charged particle in the field of an electrostatic wave packet with an arbitrary 
but finite number of harmonics. The dependence of both the edge of chaos (dissipa- 
Pu^ ' tive regime) and the deterministic diffusion on the wave packet width is predicted 



o 



theoretically and confirmed numerically, including the case of relativistic particles. 
Diverse properties of the standard maps are shown to be non-universal in the frame- 



CN ' work of the wave-particle interaction, because these maps correspond to an infinite 

> ' 

. number of waves. 



PACS numbers: 05.45.-a, 52.40.Db 



During the past quarter of century the Hamiltonian [1] and dissipative [2] versions of the 
so-called standard map (SM) have been widely studied as basic models of the dynamics of 
> ■ charged particles in the field of an electrostatic wave packet (see, e.g., refs. [3-4]): 
>< 

» I ^ ^ 

Ci ' X + 'JX = > En sin (knX — UOnt) , (1) 



^ n=-N 



where En,kn, and Un are the amplitudes, wave numbers, and frequencies, respectively, of the 
{2N + 1) plane waves, and 7, e and mg are the damping coefficient, the charge and the mass 
of the particle, respectively. Specifically, the SM describes the particular case of an infinite 
set of waves having the same amplitudes, same wave numbers, and integer frequencies. 
While its diffusion properties (7 = 0) have been shown to be non-universal in the sense 
that it does not generalize to a wave spectrum with uncorrelated phases [5], its assumption 
of an infinite and uniform amplitude distribution also seems quite restrictive since it does 
not permit one to study the sensitivity of the dynamics to changes in the wave packet 



2 



width. Physically, this sensitivity yields changes in the properties and structures of the phase 
space, as can be appreciated when comparing, e.g., the cases of two waves [6-9] and infinite 
waves [4]. In this Letter a generalized model of the wave packet structure is introduced 
to take into account such a finite-size effect on the particle dynamics. Specifically, it is 
assumed that /c„ = ko, uJn = + nAu, En = En {m) = Eq sech [nnK{l — m)/ K{m)] , where 
K{m) is the complete elliptic integral of the first kind, i.e., a sech distribution is assumed 
for the amplitudes such that the effective width is controlled by a single parameter: the 
elliptic parameter m. This specific form of E^ (m) is motivated by the following properties: 
{i) Enijn = 0) = E^SnQ, with 5nQ being the Kronecker delta, i.e., one recovers the (non- 
chaotic) limiting case of a single plane wave, {ii) En{m = 1) = Eo,\/n, i.e., one recovers 
the limiting case described by the SM. (Hi) For any m G [0, 1), one may define an effective 
number of harmonics forming the wave packet as follows. Let us choose quite freely a 
real number ( e (0, 1) such that N^ff is the largest integer satisfying E^^^^/Eq ^ then 
EjEo < C,Vn > Neff, i.e., Neff = A^e//M = [X(m) cosh"^ (1/C) /(7rX(l - m))] +1 
where the brackets stand for the integer part. Thus, for this choice of the wave packet 
structure, one has ^n=-N sin {knX — ujnt) = sin {kox — ooot) Yl!l=~N^ff En{iTi) cos (nAcut) 
and, after extending the summation from — oo to oo, eq. (1) transforms into the form 

eE 

x + ^x — sm{kQX — cuot) D{t;T,m), (2) 

rUe 

where T = 27r/Aa; is the characteristic period of the field and D{t;T,m) = 
2K{m)dn[2K{m)t/T;m\ /n [10], with dn being the Jacobian elliptic function of param- 
eter m. Note that this approximation does not imply a loss of generality at variance with 
the case of the SM. In a reference frame moving with the main wave, eq. (2) transforms into 
the equation 

— + sm^^ -Vj^ - [D{r;a,m) - l]sin^, (3) 

where ^ = koX — ujot, Qq = i^koEo/me)^^'^ , r = Qot, 5 = ^ojoUq'^, t] = ^, and a = QqT 
are all dimensionlcss variables and parameters. It is worth noting that a = Vk, where k 
is the stochasticity parameter of the SM [9]. The parameters ko,(jjQ,EQ, and Qq are held 
constant throughout. Physically, eq. (3) represents a damped pendulum subjected to a 
periodic string of finite pulses having an effective width and an amplitude controlled by m. 
Thus, model (2) permits one to study the structural stability of the system under changes in 
the width of the wave packet by solely varying the parameter m (and hence A^e//) between 
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and 1. In particular, the case of a single plane wave (m = 0) is described by a purely 
damped pendulum, while the case of an infinity of plane waves (SM: m = 1) is described by 
a delta-kicked rotator. The dependence of the chaotic dynamics' features on the effective 
number of harmonics will be considered here. 

Dissipative Regime: Consider first the case of weak dissipation (0 < 7 ^ 1). Since 
D{T;a,m) — 1 — {2 / EQ)Y^'^^^En{m) cos [2mrT/ a], eq. (3) may be regarded as a per- 
turbed pendulum (0 < 5, 77 <^ 1) for m e [0, 1) and then one can apply Melnikov's method 
(MM) [11-13,3] to obtain an analytical estimate of the edge of chaos in the parame- 
ter space. The application of MM to eq. (3) gives the Melnikov function (tq) = 
+ (167rVEo)Er=i^^^"M6„(a)sin(2n7rro/a), with = Sr] ± 2n6, 6„ (a) = 
q;~^ csch (nTT^/o;), and where the positive (negative) sign refers to the top (bottom) ho- 
moclinic orbit of the underlying conservative pendulum. It is well known that the simple 
zeros of the Melnikov function imply transversal intersections of stable and unstable man- 
ifolds (i.e., a homoclinic bifurcation occurs), giving rise to Smale horseshoes and hence 
hyperbohc invariant sets [12]. Prom (tq) one sees that a homochnic bifurcation (signi- 
fying the possibility of chaotic behavior) is guaranteed for trajectories whose initial con- 
ditions are sufficiently close to the separatrix of the underlying conservative pendulum 
if \2rj±7rS/2\ < U {m,a) = {4'K^/Eo)Y^^^-^n'^En{m)bn{oi), where U{m,a) is the chaotic 
threshold function. It is straightforward to obtain the following properties: (i) U{m,a) 
increases, as a function of m, as m is increased, i.e., the possibility of chaos increases as 
the spectral width is increased; {ii) U(m,a — > 0, 00) = 0, i.e., for any spectral width, the 
possibility of chaos diminishes when the small-amplitude frequency of the non-perturbed 
equivalent pendulum {d?^/dt^ + VL^ sin^ = 0) is much higher or much lower than the charac- 
teristic spectral frequency; {Hi) U {m, a) presents a maximum, as a function of the param- 
eter a, at a-aiax = ctmax (?7t.) , Vm e (0, 1), such that ccmax ("^) is a monotonously increasing 
function and exhibits the asymptotic behavior CKmax (?77. — > 1) ~ In (1 — m)""*^^^ ; and {iv) 
U{m — > l,a) = 47r^ (a) > U{m,a),\/m e [0,1), i.e., it is expected that the limit- 
ing case m = 1 be maximal with respect to the extension of dissipative chaos in parameter 
space. Note that properties (m) and {iv) mean that the dissipative SM is non-universal in 
the framework of the wave-particle interaction, because this map corresponds to an infinite 
set of waves having the same amplitudes, and hence it does not present the aforementioned 
maximum. Extensive Lyapunov exponent (LE) calculations of eq. (3) are coherent with 
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properties {i)-{iv). One typically finds that the maximal LE, A"^, presents a maximum as 
a function of a at = (m) and that both A"^(q;^3^) and increase with the 

spectral width in accordance with the MM-based predictions (It is worth mentioning that 
one cannot expect too good a quantitative agreement between amax(^^) and amax(^) be- 
cause MM is a perturbative method generally related to transient chaos, while LE provides 
information concerning solely steady chaos [12]). Figure 1 shows an illustrative example for 
three increasing values of the spectral width. 

Dissipationless Regime: When dissipation is negligible (7 = 0), system (3) is generated 
by the Hamiltonian H{^, p^,t) = p|/2 — D (r; a, m) cos = kox/flo, and MM provides an 
estimate of the width of the stochastic layer generated around the unperturbed scparatrix 
[14]: d{m,a) = U{m,a)/2 = {27i^ /Eo)'^'^^^n'^En{m)bn{c(). Thus, the aforementioned 
properties of the chaotic threshold function hold for the width of the stochastic layer and 
hence one expects an intensification of the deterministic diffusion as the spectral width is 
increased. Figures 2(b, d, f) provide an illustrative sequence for three increasing values of 
m at q; ~ OiuiaxiiTi), respectively. Moreover, one typically finds the gradual disappearance 
of the invariant curves region inside the separatrix cell (cf. figs. 2(c, b, d, f)). Numerical 
simulations also confirm that the stochastic layer exhibits a maximal width as a function of 
q; at a value close to Oimax {m')- This can be appreciated in the sequence of figs. 2 (a, c, e). 
Numerical simulations also indicate that another difference with respect to the SM is that 
the phase space of the Hamiltonian H{^,p^,t) is bounded by Kolmogorov- Arnold- Moser 
(KAM) tori at any values of a and m e (0, 1), i.e., global stochasticity is not possible for 
any arbitrary but finite spectral width (notice that, contrary to the SM, the phase space of 
the Hamiltonian H{C,,p^,t) is not periodic in p^). One thus concludes that, in the context 
of the wave-particle interaction, such a transition to global stochasticity is a peculiarity of 
the SM. 

Relativistic Regime: Similarly to the case of the relativistic standard map (RSM) 
[15,16], the relativistic generalization of model (2) (7 = 0) is necessary if the accelera- 
tion of particles is sufficiently large. The relativistic equations corresponding to model (2), 

X = pc^ {itlIqC^ + p^c^) p = —cEq sin {koX — Uot) D{t; T, m), where p, rrieo, and c arc the 
momentum and the rest mass of the particle, and the light velocity, respectively, may be 
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conveniently transformed into the dimensionless form 



-r, 



1 + ep| 



dr 



sin^D (r; a, m) , 



(4) 



where e = (Qo/^o) /(^ is the relativistic parameter, = kop/{meo^o) is the dimension- 
less momentum, and F = cuq/^q. To also characterize the relativistic dynamics in the 
coordinate- velocity phase space, {^^d^/dr), eq. (4) is rewritten as a second-order dif- 



fect is the reduction of the system symmetry: eq. (4) presents two mirror symmetries 
with respect to ^ = and — F, respectively, for non-relativistic particles, while it 
solely presents the former symmetry for relativistic particles when F > 0. Another rela- 
tivistic effect is the modification of the fixed points existing in the classical (Newtonian) 
regime: = {(O, F/Vl — ^r^) , (±7r, F/Vl — ^r^)}- Two limiting cases may be dis- 

tinguished: e <^ F~^ and e < F~^. Physically, the resonance condition e = F~^ means 
that the phase velocity of the main wave (uo/ko) is equal to the velocity of light. In 
the former case, when the phase velocity differs significantly from the velocity of light, 
it is found numerically that the stochastic motion is restricted to the neighborhood of 
the fixed points, such that the corresponding chaotic layers are bounded by KAM tori 
at any values of a and m e (0,1) (see figs. 3 and 4). For relativistic particles just 
beyond the Newtonian regime (e > 0), it is found numerically that the extension of the 
stochasticity regions in the coordinate-momentum phase space for relativistic and nonrel- 
ativistic particles, respectively, is approximately the same, while it diminishes relatively 
in the coordinate-velocity phase space for relativistic particles, as can be appreciated in 
the examples shown in figs. 4 and 3, respectively. Clearly, this shrinkage is a conse- 
quence of the relativistic constraint \d^/dT\ ^ e~^/^ — F (i.e., \x\ ^ c). Note that a 
proper comparison of the structures appearing in the two phase spaces would require 
taking the initial momentum as p^{0) = [((i^/(ir)(0) + F] |l — e [((i,^/(ir) (0) + F]^} 
where (d^/dr) (0) is the corresponding initial velocity. A preliminary quantitative es- 
timate of such a relativistic effect in the coordinate-velocity phase space can be ob- 
tained from the equation giving the first-order relativistic correction: d^^/dr'^ -|- sin^ = 
— [D (r; a, m) — 1] sin^-|- |eD (r; a, m) {d^/dr -\- F)^ sin^-|- O (e^). Similarly to eq. (3), this 



ferential equation: d'^^/dr'^ 



D (r; a, m) sin ^ [1 - e {di/dr + F)^] ' . A relativistic ef- 
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equation may be considered as a perturbed pendulum for m e [0, 1) and, after apply- 
ing MM to it, one straightforwardly obtains the following expression for the width of the 
stochastic layer: (m, a, e) = U{m,a,e)/2 = (27r^/£'o) ^^^^ n^£'„(m) [6,„ (a) — er„(«)], 
where r„ (a) = [2«^^ + n^7r^a^^/2] csch {mr^/a) and f/i?(m, a, e) is the first-order relativis- 
tic threshold function. One obtains that, for fixed m, a, the relativistic width du decreases 
as e is increased, and that this decrease is ever more noticeable as m is increased. Numerical 
results confirm these two predictions, as in the example shown in fig. 3. Proceeding similarly 
to obtain a first-order relativistic correction in the coordinate-momentum phase space by 
expanding the square root in eq. (4), the MM yields a null correction to the Newtonian 
width because the corresponding integrals vanish. Numerical results confirm this prediction, 
as in the example shown in fig. 4. As the phase velocity approaches the velocity of light 
(e < r~^), the momentum of the fixed points increases continuously, and the KAM torus 
limiting the stochastic region from above expands into the higher-momentum region (see fig. 
5(a)). When the resonance condition e = F"^ is met exactly, the particles can be accelerated 
to very high energies, which is indicated by the vertical traces in fig. 5(b). This mechanism 
of particle acceleration [17,18] is also observed in the RSM for the case uiq = 27il/T, with / 
being an integer [15,16]. Notice that, far from the resonance condition, the analysis of the 
RSM indicated [15,16] a relativity- induced decrease of the stochasticity in the coordinate- 
momentum phase space. However, the above findings clearly show that such a property does 
not hold for finite wave packets. 

In sum, a generalized model of the dynamics of a charged particle in the field of an elec- 
trostatic wave packet with an arbitrary but finite number of harmonics has been presented. 
The dependence of both the chaotic threshold (dissipative regime) and the deterministic 
diffusion on the spectral width was predicted theoretically and confirmed numerically, in- 
cluding the case of relativistic particles. Several properties of the SM and the RSM were 
shown to be non-universal in the framework of the wave-particle interaction, because these 
maps correspond to an infinite number of waves. 
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A. Figure Captions 

Figure 1. Maximal LE vs parameter a for three m values: 1-10~^, 1 — 10~^, and 1-10~^, 
which correspond to Nf.ff — 5, 7, and 8, respectively, for C, — 0.05. Vertical arrows indicate 
the position of the respective maxima. System parameters: 5 — r] — Q.l. 
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Figure 2. Trajectories in the stroboscopic (for t — na with n = 0, 400) Poincare section 
{Pe, vs ^/tt) for the system (3) with no dissipation [5 — r] — 0), and (a) m — 0.5, a = 1.4, 
(b) m = 0.95, a = 5.95 ~ Q;max(?7^ = 0.95), (c) m = 0.5, a = 5.3 ~ oiraa.^{fn = 0.5), (d) 
m = 1 - 10"^ a = 11.7 ~ amax(m = 1 - 10"^), (c) m = 0.5, a = 7.9, and (f) m = 1 ~ 10"^^ 
a = 23.5 ~ «i„ax(m = 1 - 10"^^). For fixed ( = 0.05, one has Neff = 2,3,7,14 for 
m = 0.5, 0.95, 1 - 10"^ 1 - 10-^^ respectively. 

Figure 3. Trajectories in the stroboscopic (for r = na with n = 0, ...,400) Poincare 
section {d^/dr vs ^/tt) for the relativistic system (see text), and (a) m — 0.6, e = 0, (b) 
m — 0.6, e = 1/9, (c) m — 0.9, e = 0, and (d) m — 0.9, e = 1/9. System parameters: 
Q; = 5.3,r = 0. 

Figure 4. Trajectories in the stroboscopic (for r — na with n — 0, ...,400) Poincare 
section (p^ vs ^/tt) for the relativistic system (see text), and (a) m = 0.6, e = 0, (b) m = 0.6, 
e = 1/9, (c) m — 0.9, e = 0, and (d) m — 0.9, e — 1/9. The initial conditions correspond to 
those used in fig. 4. System parameters: a = 5.3, F = 0. 

Figure 5. Trajectories in the stroboscopic (for r — na with n — 0, ...,400) Poincare 
section {p^ vs ^/n) for the relativistic system (see text), and (a) e = 0.8, and (b) e — 1. 
System parameters: m = 0.5, a = 5.3, F = 1. 
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